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Abstract 



We construct a Hamiltonian which in a scaling limit becomes equivalent to one that 
can be diagonalized by a Bogoliubov transformation. There may appear simultaneously 
a mean-field and a superconducting phase. They influence each other in a complicated 
way. For instance, an attractive mean field may stimulate the superconducting phase and 
a repulsive one may destroy it. 
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Introduction 

In quantum mechanics a mean field theory means that the particle density p{x) = 
ip*(x)tp(x) (in second quantization) tends to a c-number in a suitable scaling limit. 
Of course, p(x) is only an operator valued distribution and the smeared densities 
Pf = J dx p(x)f(x) are (at best) unbounded operators, so norm convergence is not possi- 
ble. The best one can hope for is strong resolvent convergence in a representation where 
the macroscopic density is built in. The BCS-theory of superconductivity is of a different 
type where pairs of creation operators with opposite momentum ip*(k) ip*{— k) (ip the 
Fourier transform and with the same provisio) tend to c-numbers. This requires different 
types of correlations and one might think that the two possibilities are mutually exclusive. 
We shall show that this is not so by constructing a pair potential where both phenom- 
ena occure simultaneously. On purpose we shall use only one type of fermions as one 
might think that the spin-up electrons have one type of correlation and the spin-down 
the other. Also the state which carries both correlations is not an artificial construction 
but it is the KMS-state of the corresponding Bogoliubov Hamiltonian. Whether the phe- 
nomenon occurs or not depends on whether the emerging two coupled "gap equations" 
have a solution or not. This will happen to be the case in certain regions of the parameter 
space (temperature, chemical potential, relative values of the two coupling constants). For 
simple forms of the potentials these regimes will be explicitly shown. Our considerations 
hold for arbitrary space dimension. 

1 Quadratic fluctuations in a KMS-state 

The solvability of the BCS-model [IJ rests upon the observation [g] that in an irreducible 
representation the space average of a quasi-local quantity is a c-number and is equal to 
its ground state expectation value. This allows one to replace the model Hamiltonian by 
an equivalent approximating one [0] . Remember that two Hamiltonians are considered to 
be equivalent when they lead to the same time evolution of the local observables Q. 

The same property holds on also in a temperature state (the KMS-state) and under 
conditions to be specified later it makes the co-existence of other types of phases possible. 

To make this apparent, consider the approximating (Bogoliubov) Hamiltonian 

H' B = J dp {u{p)a*{p)a{p) + ^A B (p) [a*{p)a*{-p) + a(-p)a(p)]} 

= J W{p)b*{p)b{p) , (1.1) 

which has been diagonalized by means of a standard Bogoliubov transformation with real 
coefficients (the irrelevant infinite constant in H' B has been omitted) 



b(p) = c(p)a(p) + s(p)a*(—p) 
a(p) = c(p)b(p) — s(p)b*(—p) 
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with 

c(p) = c(-p) s(p) = s(-p) 

c 2 (p) + s 2 (p) = l, (1.2) 
so that the following relations hold (keeping in mind that A, W, s, c will be /^-dependent) 



W(p) = yJu*(p)+A%(p) = W(-p) 

c 2 ( P ) - s 2 (p) = u(p)/W(p) , 2c(p)s(p) = A B (p)/W(p) (1.3) 

Hamiltonian (1.1) generates a well defined time evolution and a KMS-state for the b- 
operators. For the original creation and annihilation operators a, a* this gives the following 
evolution 

a(p) -> a(p) (c 2 {p)e- iW ^ t + s 2 (p)e iW ^ t ) - 2ia*{-p)c{p)s{p) sin W{p)t 
and nonvanishing termal expectations 

(a*(p)a(p')> = *(p , ' S ' i</;) 1 



1 _|_ e /3(w(p)-M) i _|_ e-UTO-ri j 
:= 5(p-p')M (1-4) 
(a(p)a(-p')} = <?(p - p'MpMp) tanh - ^ := £(p - p') [p] (1.5) 

m = {-p}, m = -[-p] 

c and s are multiplication operators and are never Hilbert-Schmidt. Thus different c and 
s lead to inequivalent representations and should be considered as different phases of the 
system. 

The expectation value of a biquadratic (in creation and annihilation operators) quan- 
tity is expressed through (1.4,5) 

(a*(q)a*(q')a(p)a(p')) = 5(q + q')8(p + p')[q}\p}- 
S(P ~ qW - q'){p}{p'} + 5(p - q')6(p' - q){p}{p'} (1.6) 

So far we have written everything in terms of the operator valued distributions a(p). 
They can be easily converted into operators in the Hilbert space generated by the KMS- 
state by smearing with suitable test functions. Thus, by smearing with e.g. 

e -«(p+ P ') 2 -«(5+9') 2 u (p) u ( 5 ) > v e L 2 (R d ) (1.7) 
one observes that in the limit k — > oo the first term in (1.6) remains finite 

< / dpdqv(p)v(q)[p][q] < oo , 
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while the two others vanish 

lim / dpdp'e- 2Kip+p ' )2 v(p)v(p'){p}{p'} = lim «T 3/2 f dpv 2 (p){p} 2 = 0. 

Since we are in the situation of Lemma 1 in 0, we have thus proved the following 
statement 

s- lim J dpdp'V(q, q' ,p,p')e~ K ( p+p ^ 2 a(p)a(p) = J dpV(q,q',p,—p)\p] (1.8) 

for kernels V such that the integrals are finite. 

With this observation in mind, a potential which acts for k —>■ oo like (1.1) might be 
written as 

V B = k 3/2 J dpdp , dqdq'a%q)a*(q')a(p)a(p')V B (q,q\p,p')e- K{p+p ' )2 - K(q+q ' )2 (1.9) 

with VB(q,q',P,p') = — Vb(q', q,P,p') etc., in order to respect the fermi nature of a's. 
This potential has the property 

||V|| < oo for k < oo 

||y|| — > oo for k — > oo 

Despite this divergence, potential (1.9) may still generate a well-defined time evolution. 
The strong resolvent convergence in (1.8) is essential, weak convergence would not be 
enough since it does not guarantee the automorphism property 

ri(ab) = r^riib) - r^ab) = r^(a)r^(6) . 

Note that the parameter k plays in this construction the role of the volume from the 
considerations in 0. 

In the mean-field regime we want an effective Hamiltonian 

H' B = J dp [u(p)a*(p)a(p) + A M (p)a*(p)a(p)] . (1.10) 

Here the KMS-state is defined for the operators a, a* themselves and one should rather 
smear by means of 

e- K{q - p)2 - Kiq '- p ' )2 v(p)v(p') (1.11) 
instead of (1.7), thus concluding that 

s-Jim J dpdqe-^ 2 a*(q)a(p)V M (q,q',P,p') =-j dp \ M j P J^_^ , (1.12) 

with e(p) = u(p) + Am(p)- Relation (1.12) then suggests another starting potential 

V M = ^ 3/2 / dpdp / dqdq / a*(q)a*(q')a(p)a(p')V M (q, q',P,p') e ~ K ^ 2 ~ ^ (1.13) 

with the same symmetry for the density Vm as in (1.9). However, in both cases a Gaussian 
form factor in the smearing functions (1.7), (1.11) has been chosen just for simplicity. In 
principle, this might be functions which have the 5-function as a limit. 
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2 The model 

Consider the following Hamiltonian 

H = H kin + V B + V M , (2.1) 

where H kin is the kinetic term and Vb,Vm are given by (1.9), (1.13). The first potential 
term describes the superconducting phase, similarly to the BCS-model, while the second 
corresponds to the mean field regime. As already mentioned, both these terms diverge 
in the limit k — > oo. The solvability of the model for k — > oo depends on whether or 
not it would be possible to replace (2.1) by an equivalent Hamiltonian that might be 
readily diagonalized. Remind that by equivalence of two Hamiltonians an equivalence of 
the time evolution of the local observables they generate should be understood. Therefore, 
the object of interest is the commutator of, say, a creation operator with the potential. 
With (1.8), (1.12) taken into account, it reads 

[a(k), V] = 2jdp {c(p)s(p) [p] V B (k, -k,p, -p)a*(-k) + V M (p, k,p, k) {p} a(k)} (2.2) 

The Bogoliubov-type Hamiltonian for our problem should be a combination of (1.1) and 
(1.10), that is of the form 

H B = J dp^a*(p)a(p)[Lj(p) + A M (p)} + ^A B (p)[a*(p)a*(-p) + a(-p)a(p)]^ (2.3) 

This Hamiltonian becomes equivalent to the model Hamiltonian (2.1), provided the com- 
mutator [a(k), H B — H kin ] equals (2.2). Thus we are led to a system of two coupled "gap 
equations" 

Ia mW = / V „(M { 1+ ;;,^,_ > „ + 1+e ;'X,. ri } * (2.4) 

A B (k) = /V.ftri^t.h g^ ^ (2.5) 

with 

W(p) = y/[u(p) + A M (p)]* + Al(p) . (2.6) 
c (and thus s, eq.(1.2)) are determined by either of the following conditions 

c\p) - s\p) = [u(p) + A M (p)]/W(p) , 2c(p)s(p) = A B (p)/W(p) . (2.7) 



The temperature and the interaction-strenght dependence of the system (2.4-7) encode 
the solvability of the model. 
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3 Solution to the coupled gap equations 

For general potential densities Vm, Vb solutions of the system (2.4-6) cannot be explicitly- 
written. In both low- and high-temperature limits a non-trivial "mean-field gap" is 
possible, while for the superconducting phase, characterized by a non-vanishing A#, a 
critical temperature exists beyond which such a phase is no longer possible. 

However, for some simple though reasonable potentials solvability of (2.4-7) and the 
behaviour of the solutions can be discussed in more detail. 

Thus, for the special case of an interaction concentrated about the Fermi surface and 
being constant therein, potential densities can be chosen as 

V B , M (k,p) = A S) M<S(k)<S(p) 

with 

5(k) = i [0(|k| - ^ + e) - 6(|k| - ^ - £)], 

where Q(x) is the Heaviside function. With the additional assumption uj{p) = p 2 , the 
system (2.4-6) transforms for e — > into 

l A (u) - A f | 3 *M ] flll 

2 M(rf " Xm { 1 + e emVP)-ri + 1 + e -^(VP)-M) / (3 ' 1} 

H7 (v ^) = A g tanh g^g^ or A, = (3.2) 



= \/\p + A Af (/i)] 2 + A|(/i) (3.3) 

with subsidiary conditions (2.7) correspondingly modified. 

We shall always assume W(^fji) > 0. This is not really a restriction, since the opposite 
situation might be similarly treated after performing the exchange b* <-> 6, also we might 
take fi + Am > 0. There are four energies Am,Ab,//, /3 = 1/T involved. The system 
exhibits severe dependence on their relative values, in particular, the following holds: 

(i) for all values of Am, As,// there is a purely mean-field solution 

o X 

A B = 0, W(^1)-» = A M = YT ^, 

but this will not be considered further; 

(ii) W(yfp) < \\b\- Furthermore, the condition W(yJ~p) > implies 

\ B > W(^) > /x, A B < W(v^I) < 

(iii) eq.(3.1) tells us that sign Am = sign Am and |Am| < 2|Am| ; 
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(iv) W > /j + Am and W > A#. In addition, eq.(2.7) brings a restriction on the mixing 
angles in the Bogoliubov transformation, if G [— 7r/4, 7r/4] U [37r/4, — 37r/4], so that 
y/2/2 < \c(p)\ < 1. 

Thus nontrivial solutions are only possible in the following regions: 

(a) \ B > : A B > // (the area A + U i3+ U C+ on Fig. 5); 

(b) X B <0: -2fi< A M < - ^j^ 
(limitations for the area U £>~ on Fig. 5). 

Therefore two general cases have to be distinguished, corresponding to attractive or 
repulsive superconducting potential, As<0orAB>0, that are qualitatively represented 
for increasing values of |A#| on Fig. 1. 




I. \ B > 

In this case necessarily only excitations with energies > /i may be present. In both 
limits W 3> /i + 2T and W <C // + 2T solutions can be obtained analytically. 

I.A < W(y/fi) - n > 2T, tanh (x) -> 1 

In this regime one finds 

W 7 = A B (3.4) 

Am = Amt — ~~r~ (3-5) 
*b + Am 

A b = ± Ag v /( A ^-^)( A ^ + / i + 2A ^) ( 3 - 6 ) 
As + Am 
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Thus a restriction has to be fulfilled 



A B + /i + 2A M >0, (3.7) 
which is always the case by Am > and also by negative Am, if in addition 

|A M | < (A B + /i)/2, 

that determines the area A + UB + on Fig. 5. Then also the mean-field gap energy becomes 
negative, though, as already mentioned, it is possible to demand positivity for the "pure" 
mean-field energy, ji + Am > 0. 



LB < W(y/p) - // < 2T, tanh (x) -> x 
Here the solution reads 



W M = ^ (3.8) 
Am = A M (3.9) 

Ar = =b 



AL>/ ' J (a* + Am) 2 (3.10) 



\ (A fl - 2T) 



and exists in the temperature interval 



XmXb <T«^^. (3.11) 



2(/i + Am) 

The above restriction substanciates the idea of this limit as being valid at sufficiently 
high, but nevertheless not at extremely high temperatures. 
With an accuracy of 10~ 5 one can then estimate 

Ab-4/x 

*M < ^ , 

so this asymptotic regime, when admissible, corresponds to the area B + U C + on Fig. 5. 

In the intermediate region an interesting phenomenon occurs that might be visualized 
by the plot of the two sides of eq.(3.2) — Fig. 2. There are three different sets of values 
of the parameters of the theory, for which one of the following situations is realised (see 
also Fig. l(i)): 

(a) the co-existence of the mean-field regime and the superconducting phase is not 
possible (A M ^ 0, A B = 0); 



(b) for fixed values of the parameters such a mixed phase is brought in to being and is 
uniquely determined (A M ^ 0, A B ^ 0); 
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Figure 2: Plot of both sides of eq.(3.2) for \ B > : (i) /i = 0.2 X B ; (ii) ft = 0.9 X B . 
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Figure 3: (i) "Thermalized" energy f3W(y/jT)/2 for the equilibrium chemical potential, 
compared to the spectra for potentials \i = k\ B , k = 0.1,0.3,0.5,0.7,0.9 ; (ii) Phase 
diagram for the ratio ^/W(-JJI). 



(c) a kind of bifurcation occurs, so that two mixed phases with different A M and A B 
become possible. 

More precisely, there exists a one-parameter family of equilibrium chemical potentials, 

p, e = — Arccoshy^B + Ab tanh (ArccoshyAB ) , (3.12) 
Q:=PQ/2, 

for which solution of (3.1-3) is uniquely determined. The corresponding energy is then 

found to be 

f3W(^fi) 



As tanh (Arccoshy Ab ) • 
z, 

For chemical potentials fi < fi e there is no solution at all, while for fi > fi e there are 
two solutions. The same is true also for the ratio n/W(y/Ji) — areas (a), (b) and (c) on 
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Fig. 3(ii) respectively. This, together with the comparison of the "equilibrium" energy 
to the spectra for different allowed values of the chemical potential — Fig. 3(i), shows 
that for a given value of the coupling constant A^, with the change of the temperature 
the system can pass from one phase to another, so that a phase transition occurs. When 
the critical value T c = As/2 is then reached, the superconducting phase (either one, 
or two co-existing such phases) is destroyed again. Furthermore, the two co-existing 
solutions correspond to two different mixing angles, so are obtained through two different, 
hence unequivalent, Bogoliubov transformations and this does not directly afflict the local 
stability of the physical system. 

II. \ B < 

According to (3.2), for negative values of the coupling constant X B the excitations are 
necessarily with energies W < fi. Thus, the mean-field gap must be negative, Am < 0, 
that presupposes solvability of the model only when the second coupling constant Am 
is also negative. As one sees on the plot in Fig. 4, there is no nontrivial phase structure 
present in this regime and the solution is always uniquely determined (see Fig. l(ii)). 




Figure 4: Plot of both sides of eq.(3.2) for X B < 0: (i) jj, < \ B ; (it) \i > \ B . 

II.A < W(y/ft) <^/i-2T 

The solution looks like the one for positive \ B , eqs. (3.4-6), only W = \\ B \ - However, 
the relation between the parameters changes 

A B + // + 2A M <0 (3.13) 

With the assumption fi + A M > this also means X B < —ft — A M but this does not 
strengthen the general restriction (b). 
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II.B < W(y/p) > n - 2T 

In this case the model is solvable, eqs. (3.8-10), within the temperature interval 

and under the same assumption as in the case Xb > 0, one concludes that 

Xb + 4/i 

Am < ^ (area o on r lg. 5) 

The above restrictions for the allowed regions in the (X B , AM)-plane, are depicted on 
Fig. 5. Nontrivial solutions are possible in the area A ± U B ± U C + . In the ./^-regions in 
the parameter space only a highly excited system exhibits a superconducting behaviour, 
while in the ^-regions this becomes possible also for a system whose energy is close to 
the chemical potential and in a restricted temperature interval. In both cases, when all 
parameters fixed, solutions are uniquely defined. 




Figure 5: Allowed regions in the (Xb, AM)-plane for given /i. 



4 Conclusion 

Our model has four parameters, Am, Xb, H,T , but by scaling only their ratios are es- 
sential. For infinite temperature (3 = eqs. (3. 1-3) admit only the mean field solution 
Ab = , Am = Am, W = /i+Am- This appears also from Fig. 5: since the phase structure 



11 



of the model and the very existence of solutions depend on the ratios /x/T, Xb/T, Xm/T, 
the area around the origin, so no superconducting phase present, reflects the situation 
for T — > oo. By lowering the temperature one meets also the BCS-type solution but in a 
rather complicated region in the 3-dimensional parameter space. 

However, the solutions themselves in the limiting cases I. A, II. A are temperature 
independent while in I.B, II. B they do depend on T. Therefore, for given values of 
Xm,Xb and \i the type- A solutions do not change upon heating whereas the type-B 
ones may appear and disappear, thus bringing the system through regions (a), (b) and 
(c) on Fig. 3(h). 

Note that the limiting time evolution depends on the state. This would not be the case 
if norm convergence would be present as shown in |J , but we have only strong resolvent 
convergence at our disposal. 

Whenever A# is positive, it must be also > \i. Also for negative Xb, Am and Am > — H 
there exists a finite gap for Xb- A perturbation theory with respect to As is in general 
doomed to failure since for no point on the Xb = axis there is a neighbourhood full of 
the A B 7^ phase. 

It is interesting that without a mean field (the Am = axis) there are superconducting 
solutions only for A^ > \i- An attractive mean field (Am < 0) stimulates superconductivity 
since then it also appears for negative Xb- However, too strong mean field attraction 
destroyes it again. 
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